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ABSTRACT 

Context. The high spot coverage of young active stars is responsible for distortions of spectral lines that hamper the detection of 
close-in planets through radial velocity methods. 

Aims. We aim to progress towards more efficient exoplanet detection around active stars by optimizing the use of Doppler imaging in 
radial velocity measurements. 

Methods. We propose a simple method to simultaneously extract a brightness map and a set of orbital parameters through a tomo¬ 
graphic inversion technique derived from classical Doppler mapping. Based on the maximum entropy principle, the underlying idea 
is to determine the set of orbital parameters that minimizes the information content of the resulting Doppler map. We carry out a set 
of numerical simulations to perform a preliminary assessment of the robustness of our method, using an actual Doppler map of the 
very active star HR 1099 to produce a realistic synthetic data set for various sets of orbital parameters of a single planet in a circular 
orbit. 

Results. Using a simulated time series of 50 line profiles affected by a peak-to-peak activity jitter of 2.5 kms -1 , in most cases we 
are able to recover the radial velocity amplitude, orbital phase, and orbital period of an artificial planet down to a radial velocity 
semi-amplitude of the order of the radial velocity scatter due to the photon noise alone (about 50 ms -1 in our case). One noticeable 
exception occurs when the planetary orbit is close to co-rotation, in which case significant biases are observed in the reconstructed 
radial velocity amplitude, while the orbital period and phase remain robustly recovered. 

Conclusions. The present method constitutes a very simple way to extract orbital parameters from heavily distorted line profiles 
of active stars, when more classical radial velocity detection methods generally fail. It is easily adaptable to most existing Doppler 
imaging codes, paving the way towards a systematic search for close-in planets orbiting young, rapidly-rotating stars. 
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1. Introduction 

The numerous detections of exoplanets reported in the last two 
decades highlight the relatively small population of planets dis¬ 
covered around young stars. These lack of detections, because 
of the limitations of existing techniques, are especially obvious 
in studies based on radial velocity (RV) analysis, with a very 
low numb er of planets reported around stars younger t han about 
300 myr ([Setiawan et al. 2007[ |Borgniet et al.||2014| . The en¬ 
hanced a ctivity of young Sun-like stars, linked to their fast rota¬ 
tion rate (Gal let & Bouvier|2013|), is the main barrier t o RV mea¬ 
surements ( |Saar & Donahue|1997HQueloz et al.|2001| ). Since the 
heavy spot coverage of these young objects generates a RV jit¬ 
ter that can hide most planetary signatures, filtering the activity 
noise is a prerequisite to progress towards a more accurate search 


for planets orbiting young stars. Our main motivation is to bet¬ 
ter constrain the formation mechanisms of planetary systems by 
comparing th e observed prevalenc e of hot jupiters versus model 
predictions ( [Mordasini et al.|2009| , as well as to get a more com¬ 
prehensive view of the evolution of accretion disks and debris 
disks in the presence of planets pda & Lin|2005[ ). 


A number of studies have been dedicated to progress towards 
more efficient filtering of activity jitter in RV curves (e.g. 
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imaging (DI) maps can be used to perform a RV filtering tak¬ 
ing a complex distribution of active regions over the stellar sur¬ 
face into account (Donati et al.] |2014[ ). We propose a different 
approach, which does not use the two-step method of previous 
studies in which Doppler mapping is performed prior to the RV 
analysis. Instead, we implement the recovery of orbital parame¬ 
ters in the DI inversion procedure itself. We therefore do not op- 
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Fig. 1 . Spot occupancy on the artificial star, in flattened polar 
view, from observations of HR 1099 secured in 1998 January. 
Concentric circles represent, starting from the centre, parallels 
of latitude +60, +30, 0, and -30°. Radial ticks around the star 
illustrate the rotational phases, with steps of 0.05 rotation cycle, 
following |Donati et ah] ( 2003b| ). 


erate an explicit filtering of the RV curve (we actually do not ex¬ 
tract any RVs from our data), but simultaneously get a brightness 
map and a set of orbital parameters in the output of a slightly 
modified DI model. Our method shares many similarities with 
a technique previously imp lemented to measure stellar rotation 
and differential rotation ( |Donati et al.||2000t |Petit et ah]|2002| ), 
which proved to behave very ro bustly in a number of pra ctical 
cases (|Donati et al.|2003a[|Barnes et al.|2005[|Petit et al.|2008[ 
|2Ql0l ). 

In the rest of the paper, we first present our set of simulated 
data, detail our maximum entropy approach to extract the orbital 
elements of an artificial planet, and then present some prelim¬ 
inary tests of the method for various sets of orbital parameters. 
We finally list our conclusions and mention future tests that need 
to be performed to fully assess the efficiency of this technique. 


sume a spectral resolving power R = 65,000 (with a Gaussian 
instrumental profile), and compute line profiles projected on ve¬ 
locity bins of 1.8 km s _1 for a total of 50 rotational phases evenly 
spread over three consecutive stellar rotation cycles (Fig. [2]). We 
finally simulate photon noise by adding a Gaussian noise to the 
line profiles, resulting in a S/N (signal-to-noise ratio) of 1,000 
per velocity bin. 

The resulting line profiles are consistent with typical obser¬ 
vations of rapidly-rotating, solar-type stars, and, in particular, 
young Suns. The S/N is at a level commonly reached in a num¬ 
ber of published DI studies (e.g. Donati et al. 2003b| ), especially 
when cross-correlation pseudo-line profiles are modelled instead 
of individual spectral lines (Donati et al. 1997] ). In this case, 
S/N= 1,000 in the cross-correlation profile implies a S/N about 
30 times lower in the initial spectrum. Gathering 50 rotational 
phases in a single run and in a time span sufficiently short to pre¬ 
vent any significant spot evolution is demanding but not unreal¬ 
istic, in particular, for stars with rotation periods of a few days, 
displaying significant phase evolution within a single observing 
night. We stress however that the regularity of the phase cov¬ 
erage proposed here is unlikely to be reached from the ground 
because of observation gaps during the daytime. The impact of 
these daytime observations gaps can be at leastpartly reduced in 
the case of multi-site campaigns (e.g. Petit et al. 2004| ). 

The rms RV jitter for this fake data set is equal to 0.8 km s _1 , 
as estimated from the standard deviation of the first moments 
of the line profiles (about 2.5 kms' 1 peak to peak, see Fig. [ 3 ]). 
We note that the Gaussian noise alone is responsible for a RV 
rms scatter of 66 ms" 1 only (this value being a function of the 
adopted S/N and velocity bin size), so that the standard deviation 
of RVs is dominated by the simulated spot activity. 

We finally alter the line profiles to incorporate RV shifts due 
to the presence of a close-in planet, assuming that the instru¬ 
ment benefits from perfect RV stability. For simplicity, we as¬ 
sume here a single planetary companion, following a circular 
orbit. The orbit and stellar equator share the same plane. The 
time-dependent RV shift induced by the companion is defined 
by the ratio of orbital over (stellar) rotation periods (P or b/Prot), 
its semi-amplitude ( K ), and its phase delay (0) compared to stel¬ 
lar rotation. We define the phase delay so that it is equal to zero 
if the first observed stellar rotational phase (also set to zero) is a 
time of maximal planetary-induced RV shift. 


2. Simulated data 

The spotted surface of the artificial star considered here is based 
on a Doppler map of the very active star HR 1099, reproduced 
on a sphere divided into a grid of 20,000 rectangular pixels of 
approximately equal area (Fig. [I). This map was initially recon¬ 
structed from data collected at the Anglo-Australian Telescope 
in 1998 January (Donati et al.' 2003bf and features a number of 
large spots distributed from the equator to the observed pole. 
With a total spot coverage of about six per cent, this specific 
brightness distribution is a good illustration of the high spot cov¬ 
erage of cool, rapidly rotating stars (e.g. Donati et al. 2003b[ 
IMarsden et al.|20Q6l|20TTl|Waite et al.|2015])~ 

We use this spot distribution to produce a set of synthetic 
cross-correlation profiles, using the set of tools previously em¬ 
ployed by Petit et aL|( |2002| ). We assume here that the projected 
rotational velocity (v sin i) of the star is equal to 40 km s -1 , and 
that its inclination angle is equal to 38° (following the values of 
Donati et al. |2003b| ). The unspotted surface temperature is set to 


5,500 K, and the spot temperature is equal to 3,500 K. We as- 


3. Simultaneous maximum entropy inversion of the 
brightness map and orbital parameters 

3.1. Method outline 

We use here a maximum entropy approach to simultaneously 
reconstruct the stellar brightness map and orbital parameters 
from the set of simulated observations. With this strategy, or¬ 
bital parameters are directly included in the DI model, similarly 
to relevant stellar parameters (e.g. stellar inclination or vsin /). 
Following the very basic principle of maximum entropy image 
reconstruction, the most likely value for the triplet of orbital pa¬ 
rameters is the value that produces a brightness map with the 
lowest information content, i.e. the reconstructed map featuring 
the lowest spot coverage. 

In practice, we adapt the method successfully implemented 
by Petit et al. ( |2002| ) to recover the surface differential rotation 
of active stars, by computing a grid of Doppler images spanning 
a range of values of K , P 0 rb/Pmu and 0. This method is directly 
inspire d fr om the ’’entropy landscap e” method of |Donati et al.| 
(2000]) and | Watson & Dhillon| ( |200l| , although transposed to a 




















































P. Petit et al.: Maximum entropy detection of planets 


3 



Fig. 2. Simulated line profiles (normalized to the continuum 
level) generated from the artificial star of Fig. [T] in the absence of 
a planetary companion. Successive profiles are vertically shifted 
for display clarity and the rotational phase of simulated obser¬ 
vation is given in the right side of the plot (with the integer part 
indicating the rotation cycle number). 
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Fig. 3. Radial velocity time series extracted from our set of sim¬ 
ulated cross-correlation profiles using the first moment method 
(top panel, orange dots), plotted together with radial velocities 
produced by the Doppler imaging code (blue dots). This model 
does not include any orbiting planet. The bottom panel shows 
the RV residuals. 


3D parameter space. We reconstruct the brightness maps using 
the DI code of Donati & Collier Ca meron ( 1997|), wh ich im¬ 
plements the two-component model of Collier Cameron] ( |1992| ). 
This model is based on the assumption that each pixel of the 
reconstructed stellar surface contains a fraction / of unspotted 
photosphere and a fraction 1 -/ of dark spot. From the same sim¬ 
ulated data set (i.e. for the same simulated value of K , P 0 rb/Prot 
and 0), we reconstruct a set of 1,000 to 3,375 Doppler maps, 
assuming different values of the ( K , Porb/Prot^) triplet in the in¬ 
version process. We derive from each reconstructed map a value 
of the fractional spot coverage and therefore obtain a 3D spot¬ 
tedness landscape in the ( K , P OY b/Prot ,0) space (Fig. [4]). To de¬ 
termine the reconstructed value of ( K , F or b/F ro t,0), we adjust a 
3D paraboloid in the ( K , P O rb/^iot>0) space, limiting the fit to a 
few tens of points around the point of minimal spottedness. By 
locating the point of minimal information content of the Doppler 
map, the paraboloid fit gives the most likely output value for ( K , 
Porb/ Protifi)- 

Similar to the method adopted by Petit et al.] (2QQ2| , error 
bars shown in Fig. mm and [7] are derived from the paraboloid fit 
by determining the projection on the three axis of the A^ 2 = 1 
surface surrounding the minimum, following the prescrip¬ 
tion of Press et al. ( |1992 ). Since the reconstructed orbital pa¬ 
rameters are computed from DI models, their estimate benefits 
from an efficient noise reduction technique (the maximum en¬ 
tropy method), which tends to decrease the uncertainties in their 
reconstructed values beyond the capabilities of a simple^ 2 min¬ 
imization approach. In this sense, usings 2 values alone to derive 
statistical uncertainties, as we do here, implies that all error bars 
discussed in the following sections are likely over-estimated. 

We do not perform an explicit filtering of activity jitter in 
line profiles, as done by Donati et al. ( |2014| ). Instead, we obtain 
the orbital parameters as part of the image reconstruction itself, 
although the reconstructed line profiles produced by the DI code, 
in principle, can be used to perform a jitter filtering of observa¬ 
tional data to double check the orbital parameters identified with 
our method (Fig [3]). We emphasize that the rough atmospheric 
model used in the DI code is the same as that used to produce 
our fake data, so that we implicitly assume here that the DI at¬ 
mospheric model is perfect. 


3.2. Test case without a planet 

The first test of our technique consisted of simulating a star with¬ 
out any planet (K , Forb/^rot,0)=(0,0,0) and checking the out¬ 
come of the Doppler inversion. The simulated set of line pro¬ 
files is shown in Fig. [2] and the resulting map (not shown here) 
is reconstructed with a reduced^ 2 close to unity (^=0.97). The 
initial spot configuration is correctly recovered, with the location 
and area of the reconstructed spots consistent with the initial spot 
configurations. Using the set of line profiles produced by the DI 
code to remove the activity-induced RV fluctuations from our 
simulated data, we end up with a rms residual jitter of 68 m s -1 , 
which is very slightly larger than the rms scatter induced by the 
simulated photon noise alone (bottom panel of Fig. [3]), showing 
that our simulated observations were adjusted down to the noise 
level, with no detectable residuals of the activity jitter in the fil¬ 
tered RV curve. 

The reconstructed orbital parameters are shown as the left¬ 
most points of Fig. [5] As expected in this test case, Forb/^rot 
and 0 turn out to be unconstrained in the entropy space. Their 
error bars are extremely large, translating the fact that the en¬ 
tropy landscape is mostly flat in these two directions. The recon¬ 
structed value of K is equal to 51 ±81 ms -1 , showing that the 
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modified DI code cannot reconcile this simulated data set with 
a RV semi-amplitude higher than about « 130 ms -1 . The rela¬ 
tively large error bar on the reconstructe d K ob tained at K m = 0, 
compared to higher K m values (see Sect. 3.3.1 ), is a consequence 
of the correlation between the different parameters of the triplet, 
as observed in Fig. [4] (so that the large error bars on the estimate 
of P orb /P rot and 0 also tend to enlarge the error bar on K). 


3.3. Orbital parameter reconstruction 

The next test consisted of running a series of simulations with 
different input values of ( K , P 0Y b/Pmt , 0) to test the robustness 
of the method for different orbital configurations. 


3.3.1. Variations of K 

We first assumed a constant simulated value of (F 0 rb/^rot)in 
(equal to 1.7), as well as a constant phase delay 0i n = 0.5. We 
then varied our simulated K m value from 0 to 200 m s -1 (Fig. [ 5 }. 
We observe that the reconstructed values of K are very close to 
the input values (to within 10 per cent) as long as K m > 20 m s -1 , 
with typical error bars smaller than about 20 ms -1 (up to 35 
ms -1 at = 60 ms -1 ). For even smaller simulated values of 
K m , the reconstructed values display much larger biases towards 
higher RV amplitudes, with error bars increasing very fast when 

approaches zero. 

The reconstructed values of Forb/^rot are also very close to 
the input values (to within 2%), as long as the input K{ n is larger 
than 20 m s -1 . The error bars on P 0 rb/Pmt progressively increase 
while K{ n decreases, going from 0.05 at = 200 ms -1 to 0.6 
at = 30 ms -1 , and remain large whenever K m < 50 ms -1 . 
The reconstruction biases also become much larger for K m < 
20 ms -1 

The behaviour of the reconstructed 0 is mostly similar to that 
of P 0Y b / P mu with more systematic biases at low values of K m . 


3.3.2. Variations of F 0 rb/^rot 

We also performed a series of simulations at fixed values of 

= 100 ms -1 , and 0i n = 0.5. Here, (F 0 rb/^rot)in was varied 
from 0.5 to 2 (Fig. [6]). The reconstructed (Forb/^rot) is always 
close to its input value, without any systematic bias. The error 
bar on the recovered value of (F 0 rb/^rot) is steadily increasing to¬ 
wards larger values of (R 0 rb/^rot)in, showing that it is preferable 
to spread the observations over several orbital periods, rather 
than getting a denser sampling of one single orbital cycle. 

The most striking feature in the reconstructed K values is a 
clear bias of up to 50% when F 0 rb ~ Pmu although the recon¬ 
structed (F 0 rb /Rrot) and 0 are much less affected by this phe¬ 
nomenon. The poor reconstruction obtained for a planet close to 
co-rotation, already reported by Donati et al. ( |2014| ), illustrates 
that the modified DI code cannot fully distinguish between a spot 
and a planetary signature as long as the rotational and orbital 
modulations are not following a different period. As illustrated 
in Fig.[8j adding a different noise pattern to our data (with a same 
resulting S/N) does not remove this bias. Apart from this zone, 
biases in the reconstructed K do not exceed 20%, and get even 
better at (F 0 rb/^rot)> 1.6. 

3.3.3. Variations of 0 

The last test consisted of varying 0i n , at fixed values of K m and 
(Forb/^rot)in (Fig. [7]). The outcome is an accurate reconstruction 


of the triplet for all values of 0i n , with limited biases for all or¬ 
bital parameters, suggesting that the 0 parameter is the less crit¬ 
ical of the triplet in this reconstruction procedure. 


4. Impact of uncertainties in reconstruction 
parameters 

In our ideal numerical simulations, setting reconstruction param¬ 
eters in the DI code is straightforward in that we know the pa¬ 
rameters used to generate the artificial star and fake time series 
of observations exactly. This question is notoriously more diffi¬ 
cult with real data, and errors in the parameters adopted in the 
DI procedure are known to generate biases in the reconstruction 
of stellar brightness maps ( [Unruh & Collier Cameron] 1995 [ ). We 
carried out a preliminary investigation of the outcome of wrong 
values of DI input parameters in the reconstructed triplet of or¬ 
bital parameters. We based this ne w serie s of tests on simulations 
similar to those discussed in Sect. |3.3.21 i.e. where the RV am¬ 
plitude K is fixed and equal to 100 ms , as well as 0, which is 
set to 0.5, while F 0 rb/^rot is varied between 0.5 and 2.0. The re¬ 
constructed K values obtained in this new series of simulations 
are shown in Fig. [8] 

A first test consisted in varying the final x Y values of the 
models. Instead of the default target value of 0.97, which is set 
to be the highest^ value below which the rms of the RV residu¬ 
als cannot be reduced further,^ = 1.0 and 0.94 are successively 
adopted (circles in Fig. [8]). We observe that reconstruction biases 
in K values tend to decrease for decreasing^ values and the as¬ 
sociated error bars. However, the observed changes remain lim¬ 
ited compared to other effects. In particular, adopting^ = 0.97 
again with a different noise pattern in our data (simply obtained 
with a different seed to generate random numbers, red crosses 
in Fig. [8]) roughly leads to the same amplitude of changes in the 
biases and error bars. Similar conclusions are reached for the re¬ 
construction of P orb/ P rot and 0 (not shown here). Below x Y = 
0.94, most of our models did not converge at all. 

As a second test, we modified the inclination angle assumed 
during the inversion process, adopting i = 33 and 43° instead 
of i = 38° (triangles in Fig. [8]). While some of these models 
still provided us with reconstructed K values close to the input 
value, the biases become much larger in a number of cases to 
the point where several points fall outside of the boundaries of 
our figure. The obvious inaccuracy in the inclination angle is 
also highlighted by the best achievable x Y * which is more than 
doubled that of our models assuming an optimal value for the 
stellar inclination. In practice, these large x Y values suggest that 
offsets on the inclination angle in the DI input parameters should 
remain smaller in real cases. Not surprisingly, statistical error 
bars are also very large, up to ten times wider than those of our 
optimal model. 

Finally, we ran a first estimate of the impact of stellar differ¬ 
ential rotation on the reconstructed orbital parameters (Fig. [9]). 
We assumed here that the stellar surface is altered by a latitudi¬ 
nal shear that follows a simplified law: Q(/) = Q e q - sin 2 (/). JD, 
where D(/) is the rotation rate at latitude /, Q eq is the rotation rate 
of the equator, and d£l is the difference of rotation rate between 
the pole and the equator. We varied the dimensionless parameter 
y = d£l/£l Qq from zero (solid-body rotation) to 0.02, assuming 
(K , P orb/ P rot? 0) = (100 m s -1 ,1.7,0.5). Assuming a star rotating 
in three days, y = 0.02 roughly corresponds to 80% of the solar 
shear level. In this case, our artificial data (assumed to be gath¬ 
ered over three stellar rotation cycles) are collected over nine 
consecutive days. Some of our DI models compensate for differ- 
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ential rotation (Petit et al .|2002 orange dots). Other DI models 
simply ignore the shear by assuming solid-body rotation (green 
dots), showing the consequence of obtaining an incorrect esti¬ 
mate of the shear, or of using an incorrect type of differential 
rotation law. Fig. [9] shows that the biases and error bars remain 
roughly constant for DI models incorporating a differential rota¬ 
tion law. Up to y « 0.012, the situation is mostly the same for 
models assuming solid-body rotation. However, biases steadily 
increase past this point whenever the surface shear is ignored. 
This evolution is coming along with a ^ 25% increase of the 
Xr, and a significant increase in the statistical error bars. Above 
y = 0.02, the majority of our models based on solid-body ro¬ 
tation do not provide any well-defined spottedness minimum in 
the 3D entropy landscape, so that no reliable values of orbital 
parameters could be derived. The two differential rotation pa¬ 
rameters can be added to the orbital parameter triplet to perform 
a simultaneous optimization of all five parameters. In this case, 
the additional dimensions to be probed would likely require us to 
adopt an MCMC approach to identify the best model, instead of 
the time-consuming systematic paving of a very large 5D grid. 
Inclusion of the stellar inclination and v sin i is possible as well, 
leading to a 7D space to be explored. 


5. Conclusions 

In spite of a simulated set of line profiles displaying an activity 
jitter as high as 2.5 km s -1 peak to peak, we are able to robustly 
recover orbital elements of a fake planet with RV signatures of 
a few tens of ms -1 , i.e. of the order of magnitude of the RV 
scatter due to the photon noise of our simulations. This is true as 
long as the planet is not in a co-rotating configuration, in which 
case offsets of up to 50% are observed in the reconstruction of 
K (assuming K m = 100ms -1 ), while (F or b/Rrot) and 0 are kept 
closer to their input values. 

Our preliminary tests support the capabilities of this max¬ 
imum entropy method in a systematic search for hot jupiters 
around extremely active stars, as long as the collected time se¬ 
ries is well adapted for DI, i.e. with a dense rotational sampling 
collected over a few rotation periods. However, we need to com¬ 
plement the present set of simulations, with a more comprehen¬ 
sive set of models to determine the detailed effect of a number of 
parameters (e.g. v sin i, stellar inclination, orbital and rotational 
phase sampling, surface differential rotation) on the reconstruc¬ 
tion of planetary orbits. Varying the inclination will mostly be¬ 
come critical at low i values (i.e. large K / sin(z) values), since the 
detectability of low-mass planets will be much reduced in this 
geometrical configuration. The impact of v sin i is more complex 
to address. The present approach makes use of the fact that, at 
large projected rotational velocities, the spectral signatures of 
individual spots are bumps restricted to a fraction of the line 
profile, while the planetary signature affects the line profile as a 
whole. Low v sin i values will make this disentangling less easy; 
the limit of efficiency of the technique is therefore similar to the 
v sin i limit of accurate DI. In any case this threshold is below 
vsin i = 10 kms -1 , since Donati et al., in press, demonstrated 
that the method is still operative for v819 Tau with v sin i = 9.5 
km s -1 . On the other hand, if higher v sin i values improve the DI 
efficiency, at the same time it would reduce the RV accuracy, so 
that an upper limit in projected rotational velocity should exist 
as well. 

In addition, the optimization of our detection threshold will 
likely require the DI reconstruction of both cool spots and facu- 
lae ( fDonati et al.|2014| . Future tests will also assess the benefit 
of obtaining observations at different epochs to take advantage of 


the changing spot occupancy to improve the accuracy of planet 
detection and characterization in repeated measurements. The 
impact of non-rotationally modulated events that a classical DI 
code is not able to model, such as flares or limited spot lifetime, 
will also deserve special attention. 

We also stress that the tests performed here were limited to a 
single planet following a circular orbit. More complex configu¬ 
rations involving an orbital eccentricity or multiple companions 
will necessitate probing additional dimensions in a multi-D en¬ 
tropy landscape. Using any real data set, the average RV of the 
star-planet system is another parameter that cannot be ignored in 
the multi-D landscape (see Donati et al., in press). 

The present maximum entropy method constitutes a very 
simple way to extract orbital parameters from heavily distorted 
line profiles of active stars, paving the way towards a systematic 
search for close-in planets around young, rapidly-rotating stars. 
It is easily adaptable to most existing DI codes, and may be em¬ 
ployed to revisit a number of existing data sets initially collected 
for DI purpose alone. One specificity of this technique is the si¬ 
multaneous extraction of a DI map and of orbital elements, in a 
one-step procedure much simpler than the sequential approach 
generally proposed in this context. In addition to its simpler ap¬ 
proach, the lower amount of data manipulation involved here is 
likely to make the present strategy less prone to information re¬ 
moval. 

Acknowledgements. We are grateful to an anonymous referee for the suggestion 
of a number of new tests, which helped to clarify the assets and limitations of 
our method further. 
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Fig. 4. Top panel: example 3D entropy (spottedness) landscape, 
obtained from 1,000 DI models with different values of K , 
Porb/^rob and 0. The simulated parameter values were K = 
0.1 kms' 1 Porb/Prot = 1.7 and 0 = 0.5. Middle and bottom 
panels: example 2D slices in two orthogonal planes, extracted 
from the 3D orbital parameter space. Each slice is obtained from 
a set of 2,500 Doppler images. In all three panels, the colour 
scale shows the fractional spot coverage of the resulting Doppler 
maps, and the location of minimal spottedness shows the recon¬ 
structed set of orbital parameters. 
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Fig. 5. Simulations for a radial velocity amplitude K varying Fig. 6. Same as Fig.jiJ except for an orbital period varying from 
from 0 to 200 ms -1 . The orbital period is fixed and equal to 0.5 to 2 times the stellar rotation period. The radial velocity am- 
1.7 times the rotational period, as well as 0, which is set to 0.5. plitude K is fixed and equal to 100 ms -1 and 0, which is set to 

The dashed lines indicates the input values, while circles show 0.5. 
the reconstructed parameters. The l<x error bars derived on in¬ 
dividual parameters from the paraboloid fit are plotted as grey 
areas. 
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Fig. 7. Same as Fig. [i] except for an orbital phase varying from 
0 to 0.9. The radial velocity amplitude K is fixed and equal to 
100 ms -1 and P 0 rb/Prou which is set to 1.7. 



Fig. 8. Same as the top panel of Fig. [b] with errors introduced 
in the reconstruction parameters. Orange circles again show the 
points of Fig. [6j obtained with^ = 0.97 and i = 38°, i.e. our 
optimal model. Green and pale blue circles are obtained by forc¬ 
ing^ = 1.0 and 0.94, respectively, during reconstruction. Violet 
(resp. green) triangles show the effect of forcing i = 33° (resp. 
i = 43°). Crosses are obtained by adding another noise pattern 
to our fake data set (using a different seed to generate random 
numbers, still resulting in S/N^ 1,000). Typical error bars are in¬ 
dicated in the top right side of the plot, except for violet and 
green triangles (affected by much larger statistical error bars that 
would hamper the plot readability). 



Fig. 9. Reconstructed K values obtained with stellar models, in¬ 
cluding various levels of solar-like surface differential rotation. 
The dimensionless surface shear is equal to <7Q/Q e q, where Q e q 
is the rotation rate of the equator and Q is the difference of ro¬ 
tation rate between the pole and the equator. Orange dots and 
their error bars (grey area) corre spond to DI models compensat¬ 
ing for the differential rotation ( Petit et al. |2002| . Green dots are 
obtained with DI models assuming solid-body rotation. 






























